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ABSTRACT 

A  numerical  and  experimental  investigation  of  the  interaction  of  a 
pair  of  vortices  with  a  free  surface  has  been  undertaken.  The  analysis 
is  based  on  the  vortex-sheet  representation  of  the  free  surface  and  the 
use  of  the  appropriate  boundary  conditions.  The  experiments  were 
performed  in  a  large  basin  and  the  vortices  were  generated  through 
the  use  of  a  special  nozzle.  The  rise  of  the  resulting  Kelvin  oval,  the 
trajectories  of  the  vortices,  and  the  instantaneous  shape  of  the  free 
surface  were  recorded  on  a  video  tape  and  then  carefully  analyzed 
through  the  use  of  a  Motion  Analysis  system.  The  results  have  shown 
that  the  rise  of  the  vortices  not  only  gives  rise  to  two  scars,  with  a 
pronounced  hump  in  the  middle,  but  also,  and  more  importantly,  to  a 
three-dimensional  instability  heretofore  unknown.  The  measured  and 
calculated  vortex  trajectories  and  the  free-surface  shapes  at  the  corre- 
sponding times  and  Froude  numbers  are  found  to  be  in  reasonable 
agreement.  The  new  instability  will  form  the  basis  of  future 
investigations. 
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I.  INTRODUCTION 

This  is  part  of  continuing  basic  research  toward  the  understand- 
ing of  the  fundamental  mechanisms  and  physical  processes  underlying 
two-  and  three-dimensional  vortex/free-surface  interactions  in  homo- 
geneous, stratified,  and  sheared  media,  taking  into  account  ambient 
turbulence,  viscous  effects,  and  various  large-scale  instabilities  (sinus- 
oidal instability  and  vortex  breakdown)  for  ship-  and  submarine- 
related  hydrodynamics  in  a  real  ocean  environment. 

The  name  "remote  sensing"  was  first  used  in  the  early  1960s, 
when  it  was  recognized  that  the  term  "aerial  photography"  was  no 
longer  suitable  to  describe  new  kinds  of  images  formed  using  energy 
outside  the  visible  spectrum.  The  developments  in  1970s  (largely  by 
NASA)  gave  rise  to  the  development  of  a  serious  scientific  interest  in 
remote  sensing.  This  interest  led  to  the  routine  availability  of  multi- 
spectral  data  of  large  portions  of  the  earth's  surface.  As  far  as  the 
oceans  are  concerned,  the  interpretation  and  use  of  the  multi-spectral 
data  required  a  clear  understanding  of  the  physics  of  various  phenom- 
ena leading  to  surface  signatures  resulting  from  the  body-wake  free- 
surface  interaction  in  two-  and  three-dimensional,  sheared,  and/or 
stratified  media,  with  or  without  ambient  turbulence. 

As  far  as  the  ship  wakes  are  concerned,  the  most  significant  wake 
signature  has  been  thought,  for  nearly  a  century,  to  be  the  Kelvin  wave 


pattern,  with  its  38.9°  wake  angle  (see  Figure  1).  Observations  made 
with  remote  sensing  devices  (e.g.,  synthetic  aperture  radar)  have 
changed  this  perception  dramatically  and  have  shown  clearly  that  a 
narrow  strip  of  disturbed  wake  (scars  and  striations)  exists  for  several 
hours  behind  the  ship  (see  Figures  2a  and  2b)  with  little  or  no  sign  of 
the  Kelvin  wake.  These  features  have  been  termed  turbulent  wakes, 
although  it  is  not  clear  that  turbulence  alone  is  responsible  for  their 
appearance.  While  the  understanding  of  the  physics  of  this  phe- 
nomenon, through  experimentation  and  numerical  simulation,  has 
attracted  a  great  deal  of  attention,  a  satisfactory  model  describing  the 
essential  features  of  the  radar  images  has  not  yet  been  presented.  It 
appears  that  many  competing  and  interacting  processes  are  responsi- 
ble for  the  observed  scar  and  striation  pattern.  It  has  been  suggested 
that  vortices  produced  by  the  ship's  hull  may  be  responsible  for  the 
suppression  of  surface  waves  near  the  ship  track  and  an  enhancement 
of  the  waves  near  the  edges  of  the  smoothed  area.  Not  only  the  free- 
surface  signatures  generated  by  the  vortical  flow  but  also  the  natural 
waves  radiate  energy  away  from  the  wedge-shaped  wake  region  and 
thereby  affect  the  flow  itself.  It  appears  that  the  dark,  narrow  images 
are  visible  due  to  a  decrease  in  radar  return  (negative  spectral  pertur- 
bation) aft  of  a  ship  and  not  to  an  increase  in  radar  backscatter  along 
the  edges  of  the  wake. 

The  three-dimensional  temporal  development  of  the  surface  sig- 
natures resulting  from  the  interaction  of  an  ascending  inclined  vortex 


Figure  1 .  Kelvin  Waves  and  White-Water  Wake  of  a  Ship 


500  m 


Figure  2a.  SAR  Image  of  the  Narrow  Turbulent  Wake  of  a  Ship 


100  m 


Figure  2b.  Close-Up  View  of  the  SAR  Image  of  the  Narrow 
Turbulent  Wake  of  a  Ship 


pair  with  the  free  surface  was  first  described  by  Sarpkaya  [Ref.  1]  and 
Sarpkaya  and  Henderson  [Ref.  2].  When  a  trailing  vortex  pair  generated 
by  a  lifting  surface  approaches  the  free  surface,  with  or  without  mutual 
induction  instability  and/or  vortex  breakdown,  the  vortices  and/or  the 
crude  vortex  rings  give  rise  to  surface  disturbances  (scars  and  stria- 
tions,  see  Figures  3a  and  3b).  The  striations  are  essentially  three- 
dimensional  free- surface  disturbances,  normal  to  the  direction  of  the 
motion  of  the  lifting  surface,  and  come  into  existence  when  the  vortex 
pair  is  at  a  distance  equal  to  about  one  initial  vortex  separation  from 
the  free  surface.  The  scars  are  small,  free-surface  depressions  on 
either  side  of  a  nearly  symmetrical  hump,  directly  above  the  plane  of 
symmetry  of  the  vortex  pair.  They  come  into  existence  towards  the 
end  of  the  pure  striation  phase  and  when  the  vortices  are  at  a  distance 
equal  to  about  sixty  percent  of  the  initial  vortex  separation  (b0)  from 
the  free  surface.  The  scars  are  created  and  driven  by  the  vortices  and 
both  move  outward.  Sarpkaya  [Ref.  1]  and  Sarpkaya  and  Henderson 
[Ref.  2]  found  that  the  vortices  diffuse  rapidly  and  their  circulation 
decreases  slowly  at  first  (prior  to  the  occurrence  of  scars)  and  rapidly 
thereafter  due  to  the  counter-sign  vorticity  generated  in  the  overlap- 
ping regions  of  the  vortex  pair  in  the  recirculation  cell  and  at  the  free 
surface. 

Previous  studies  of  a  vortex  (or  vortex  pair)  in  ground-  or  free- 
surface  effect  have  been  generally  concerned  with  the  boundary  layer 
development,    generation    of   oppositely-signed    vorticity,    the    flow 


Figure  3a.  Scars  and  Striations  Resulting  From 
a  Pair  of  Trailing  Vortices 


Figure  3b.  Scars  and  Vortical  Structures  Resulting  From 
a  Pair  of  Trailing  Vortices 
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resulting  from  a  single  vortex  held  fixed  in  a  uniform  flow,  and  the  so- 
called  "rebounding"  phenomenon,  according  to  which  the  vortices 
move  at  first  toward  the  boundary,  as  on  an  inviscid  trajectory,  and 
then  away  from  it  [Refs.  3-6].  The  boundary  surface  may  be  either  a 
rigid  boundary  at  which  the  no-slip  condition  is  satisfied  or  a  clean 
free  surface  corresponding  to  zero  shear  stress.  In  either  case,  the 
boundary  surface  gives  rise  to  oppositely-signed  vorticity  and  the  net 
circulation  of  the  vortices  decreases  with  time  [Ref.  5]. 

Sarpkaya  and  Henderson's  [Refs.  1,  2]  theoretical  analysis  of  the 
scars  created  by  the  trailing  vortices  was  based  on  the  classical  solu- 
tion of  Lamb  [Ref.  7],  assuming  the  vortices  to  be  two-dimensional  and 
the  free  surface  to  be  a  rigid  plane.  For  small  Froude  numbers  F 
(F  =  V0/(/gbo,  where  V0  is  the  initial  mutual  induction  velocity  of  the 
vortices  and  b0  is  the  initial  vortex  separation),  each  scar  front  was 
shown  to  coincide  with  the  stagnation  point  on  the  Kelvin  oval,  formed 
by  one  of  the  pair  of  the  trailing  vortices  and  its  image.  For  Froude 
numbers  larger  than  about  0.15,  not  only  the  deformation  of  the  free 
surface  but  also  the  nonlinear  interaction  between  the  said  deforma- 
tion and  the  motion  of  the  vortices  are  significant.  Thus,  the  free  sur- 
face may  no  longer  be  assumed  rigid.  Even  though  it  was  fully  realized 
at  the  outset  that  the  problem  to  be  ultimately  solved  is  the  under- 
standing of  the  effects  of  three -dimensionality,  shear,  stratification, 
and  the  parameters  characterizing  the  body-shape  and  motion  on  the 
evolution  of  the  surface  signatures,  the  relative  ease  of  the  study  of  the 
two-dimensional     flow    attracted     the     immediate     attention     of 
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experimentalists  and  numerical  analysists  alike  in  the  past  few  years 
(Sarpkaya,  et  al.  [Ref.  8];  Willmarth,  et  al.  [Ref.  91;  Dahm,  et  al.  [Ref. 
10];  Marcus  [Ref.  11];  and  Telste  [Ref.  12]). 

Sarpkaya,  et  al.'s  two-dimensional  experiments  were  conducted  in 
a  12-foot-long,  three-foot-wide,  and  four-foot-deep  (12'  x  3'  x  4')  water 
basin.  The  nearly  two-dimensional  vortex  pair  was  generated  through 
the  use  of  two  counter-rotating  plates.  Neutrally  buoyant  fluorescent 
dye  was  used  to  visualize  the  flow.  Sarpkaya,  et  al.  [Ref.  8]  also  carried 
out  a  numerical  analysis  of  the  scars  through  the  use  of  vortex  dynam- 
ics. Their  results  have  shown  that  the  two-dimensional  evolution  of  the 
free  surface  can  be  calculated  (up  to  certain  normalized  times) 
through  the  use  of  line  vortices  or  dipoles.  The  analysis  becomes  more 
robust  as  the  Froude  number  is  increased.  At  larger  times,  the  results 
become  unrealistic  because  the  laminar  and  turbulent  diffusion  of  vor- 
ticity  are  ignored.  No  wave  train  was  observed  on  either  side  of  the 
scars  for  any  of  the  Froude  numbers  encountered  in  the  analysis  or 
experiments. 

Willmarth,  et  al.  [Ref.  9]  and  Dahm,  et  al.  [Ref.  10],  at  the  Univer- 
sity of  Michigan,  adopted  Sarpkaya,  et  al.'s  [Ref.  8]  idea  of  counter- 
rotating  plates  to  generate  vortices  and  carried  out  similar 
experiments  on  vortex-pair  free-surface  interaction.  They  too  have 
used  vortex  methods  to  calculate  the  path  of  the  vortices  and  the 
deformation  of  the  free  surface.  However,  no  detailed  comparisons 
were  presented  between  the  observed  and  calculated  scar 
characteristics. 
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Marcus  [Ref.  11],  working  with  Berger  at  UC  Berkeley,  used  the 
finite-difference  method  to  investigate  the  inviscid,  two-dimensional 
interaction  between  a  pair  of  counter-rotating  line  vortices  and  a  free 
surface.  He  has  not  conducted  any  experiments  and  encountered  con- 
siderable numerical-stability  problems  at  lower  Froude  numbers.  At 
higher  Froude  numbers,  his  predictions  were  somewhat  larger  than 
those  measured  and  calculated  by  Sarpkaya,  et  al.  [Ref.  8].  Telste 
[Ref.  12]  used  vortex  dynamics  to  investigate  the  same  phenomenon 
(inviscid  two-dimensional  interaction  between  a  two-dimensional  vor- 
tex pair  and  the  free  surface).  At  medium  and  lower  Froude  numbers, 
his  calculations  required  extensive  numerical  filtering  to  avoid  or  to 
delay  numerical  instability.  The  calculations  are  generally  terminated 
when  the  free  surface  forms  a  sharp  corner  or  when  difficulties  are 
encountered  in  obtaining  rapid  convergence.  There  does  not  seem  to 
be  any  direct  relation  between  the  breakdown  of  the  numerical  model, 
breakdown  of  the  waves,  and  the  attainment  of  maximum  scar  condi- 
tions. Telste  offered  no  comparisons  with  experiments. 

It  is  clear  from  the  foregoing  that  the  investigation  first  initiated 
by  Sarpkaya  [Ref.  1]  at  the  Naval  Postgraduate  School  quickly  mush- 
roomed into  a  full-scale  scientific  effort  on  vortex/free-surface 
interactions,  partly  because  of  its  intrinsic  interest  and  partly  because 
of  its  far-reaching  practical  applications  in  ship-  and  submarine- 
related  hydrodynamics  in  a  real  ocean  environment.  It  is  also  clear 
that   much    remains    to   be    accomplished,    both    theoretically   and 
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experimentally,  toward  the  understanding  of  the  evolution  of  surface 
signatures. 

The  present  investigation  has  a  two-fold  purpose: 

1.  To  conduct  further  parametric  studies  of  the  variables  affecting 
the  performance  of  the  numerical  model  through  the  use  of  the 
vortex-sheet  representation  of  the  free  surface;  and 

2.  To  investigate  in  greater  detail  the  motion  of  a  vortex  pair 
approaching  and  intruding  through  the  free  surface  and  to  delin- 
eate the  characteristics  of  the  scars  and  striations. 
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II.  NUMERICAL  SIMULATION 

A.      DESCRIPTION  OF  THE  MODEL 

The  flow  model  consists  of  two  counter-rotating  vortices  (CRVs) 
of  constant  circulation,  moving  towards  the  free  surface,  and  a  vortex 
sheet,  represented  by  a  number  of  equi-spaced  line  vortices  at  the  free 
surface.  The  model  is  a  generalized  boundary  integral /vortex  method 
developed  by  Baker,  Meiron  and  Orszag  [Ref.  13]. 

The  use  of  Bernoulli's  equation  on  either  side  of  the  fluid  interface 
leads  to 

^  +  i(u211  +  u2n,)+gn1  +  ff  =  Fl(t)  (la) 


and 


<^2         !  P2 

-5t"  +    *  (U*  +  Un2>   +  ^  2  +    p-  =  F2^  I  lb) 


in  which  <J)  is  the  velocity  potential,  u  is  velocity,  g  is  gravitational 
acceleration,  t  is  time,  r\  is  surface  elevation  above  mean  water  level,  p 
is  density,  p  is  pressure,  s  is  distance  along  surface,  the  suffixes  t  and 
n  denote  components  tangential  and  normal  to  the  surface,  and  the 
suffixes  1  and  2  denote  the  lower  and  upper  surface.  Furthermore,  on 
the  surface  one  has 

Pi=P2  =  P 
Unl  =  Un2  =  un 
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111=112  =T1 
Yt  =  Utl-Ut2 

ut  =  (uti +ut2)/2  (2) 

U  =  (ui  +  U2)/2 
V  =  (Vl  +v2)/2 
r  =yt  ds  =  ydx 

The  use  of  the  foregoing  plus  the  free  surface  condition  given  by 


^-£(*ax+<D2y)+F-2Ti  =  0 


(3) 


lead  [Ref.    14]  to  the  momentum  equation  for  the  calculation  of  the 
evolution  of  vorticity 


^  =  -^(Uy)+2At 


au    TTau    fdv    TTav^an 


>2,y 


C  y[  dy 
~5x~ 


-CSy 


^2      \ 


an 

+  g-r- 
6dx 


(4) 


where  the  Atwood  number  At  =  (p2  -  Pi)/(p2  +  Pi)-  For  the  air/water 
interface,  At  =  -1. 

The  velocity  of  a  vortex  may  be  decomposed  into  two  parts:  One 
stemming  from  the  use  of  line  vortices  (i.e.,  the  two  CRVs  and  the 
surface  vortices,  equi-spaced  at  mid-points  of  the  segments)  and  the 
other  from  the  corrections  to  the  first  due  to  the  distributed  nature  of 
vorticity  along  the  sheet.  Thus  one  has: 


Uk/V0  =  fUi  +  U2)/V0 


(5a) 
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Vk/V0  =  (Vi  +  V2)/V0  (5b) 

where  the  Ui  and  Vi  components  of  velocity  are  obtained  from  the 
complex  function  w(z).  Normalizing  the  strengths  of  all  vortices  by  the 
strength  of  the  counter-rotating  vortex  pair,  all  velocities  by 
Vo  (=  ro/27ibo),  and  all  distances  by  b0,  one  has 

w  =  Mn(z  -zj  -  ^ln(z  +  zj  -  £-^rkln(z  -zk)  +  XIrkln(z  +  zk)    (6) 
1  1  k=i l  k=i l 

in  which  z0  and  -z0  are  the  positions  of  the  CRVs  and  Zk  and  -Zk 
denote  the  positions  of  the  free  surface  vortices.  The  velocities  Ui  and 
Vi  are  then  given  by 

(Ui  -iVi)/V0  =  dw/dz  (7) 

and  the  velocity  components  U2  and  V2  are  given  by  [Ref.  14] 


U    =  ^ 


dy_  dn      Xo^  (8a) 

{dx'dx  +  2  dx2  J 

k 


k 


From  a  mathematical  point  of  view,  one  would  like  to  have  the 
free  surface  extend  from  ^0  to  +00  and  the  CRVs  originate  at  (0.5,  -00) 
and  at  (-0.5,  -«»).  This  is  impossible  from  a  numerical  point  of  view.  It 
is  necessary  to  obtain  guidance  from  experiments  to  limit  the  com- 
puter time.  Observations  have  shown  that  the  free  surface  rises  in  a 
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relatively  small  region  directly  above  the  CRVs  and  no  waves  are  gen- 
erated beyond  the  scars  (see  Figure  4).  In  fact,  the  linearized  analysis 
of  Kochin,  et  al.  [Ref.  15]  of  the  motion  of  a  line  vortex  under  a  free 
surface  has  shown  that  the  resulting  steady  wave  train  far  from  the 
vortices  would  be  negligibly  small.  The  sample  calculations  led  to  the 
conclusion  that  the  free  surface  can  be  restricted  to  a  region  extend- 
ing from  x  =  -10  to  x  =  +10.  Furthermore,  the  CRVs  are  assumed  to 
originate  at  a  depth  of  y  =  -5,  in  view  of  the  fact  that  the  free  surface 
does  not  begin  to  deform  until  the  CRVs  arrive  at  about  y  =  -1. 

Theoretically,  one  can  calculate  the  potential  function  §  from  the 
complex  velocity  potential  w,  the  velocity  of  the  vortices  from  dw/dz  = 
u-iv  (excluding  the  vortex  whose  velocity  is  being  calculated),  and  the 
elevation  of  the  free  surface  from  equations  3  or  4  for  a  given  Froude 
number,  and  trace  the  evolution  of  the  free  surface  as  a  function  of 
time.  This  relatively  simple-sounding  procedure  is  anything  but  sim- 
ple, primarily  because  the  numerical  short-wavelength  instabilities 
eventually  lead  to  large-scale  instabilities.  The  critical  time  for  the 
onset  of  these  instabilities  decreases  with  decreasing  Froude  number. 
It  is  because  of  this  reason  that  the  preliminary  calculations  were  ini- 
tially restricted  to  large  Froude  numbers.  To  be  sure,  the  instability  is 
not  limited  to  boundary-integral  methods.  It  is  also  an  inherent  char- 
acteristic of  the  space  discretization  methods.  In  fact,  the  choice  of  a 
particular  method  depends,  more  or  less,  on  the  user's  past  experi- 
ence,   the    availability   of  a   computer   code,    the    objectives    of  the 
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Figure  4.  Photograph  of  a  Typical  Free-Surface  Deformation 
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investigation,  and  the  particular  problem.  A  lucid  discussion  of  the 
advantages  and  disadvantages  of  various  techniques  is  presented  by 
Yeung  [Ref.  16].  Previous  investigators  dealing  with  marching  methods 
for  nonlinear  water  wave  problems  have  used  various  filtering  [Ref.  17] 
or  regridding  techniques  to  eliminate  or  delay  the  instability.  In  gen- 
eral, it  is  not  easy  to  assess  the  accuracy  of  the  results.  Sensitivity 
analysis,  parametric  studies,  and  comparison  with  experiments  may 
help  to  assess  the  error  introduced  by  numerical  smoothing  and  fil- 
tering techniques  [Ref.  18].  One  must  also  bear  in  mind  that  the 
numerical  model  might  not  necessarily  correspond  to  a  physically 
realizable  flow  (turbulent  diffusion,  axial  instabilities,  wall  effects,  etc.). 
In  the  present  investigation,  it  was  decided  to  use  Rosenhead's 
[Ref.  19]  smoothing  technique,  according  to  which  any  complex  veloc- 
ity is,  in  general,  desingularized  by  introducing  a  cut-off  distance  5, 
given  by  (see  also  Sarpkaya  [Ref.  18]): 


1   v*   r' 

U    -IV,    =  77: —  >     - — 
J         J       2irc  /   jZi  - 


(2 

(9) 


Z.        ,  _      ,2     .     c2 

k      z 


;j-zkf +5' 


The  use  of  other  smoothing  techniques  (e.g.,  the  Longuet-Higgins 
and  Cokelet  technique  [Ref.  17])  was  rejected  because  it  tends  to  over- 
smoothen  the  free  surface.  Furthermore,  as  argued  by  Roberts 
[Ref.  20],  the  five-point  scheme  of  Longuet-Higgins  and  Cokelet  may 
lead  to  hopelessly  large  errors  in  some  applications.  According  to 
equation  6,  the  solution  should  approach  the  inviscid  flow  solution  as 
5  — >  0.  Numerous  calculations  have  shown  that  a  8/ Ax  value  between  1.1 
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and    1.5  yields  nearly  identical  results.  The  final  calculations  were 
performed  using  8/ Ax  =  1.25. 

The  specific  details  of  the  numerical  calculations  are  as  follows:  (i) 
calculate  the  strength  of  the  vortex  sheet  and  the  strength  of  the  sur- 
face vortices  by  iteration  until  3-y/3t  and/or  d\J/dt  and  dV/dt  converge  to 
an  acceptable  limit,  (ii)  find  the  velocity  of  all  vortices  (including  those 
of  CRVs)  and  convect  them  for  a  time  interval  At  using  a  second  order 
modified  Eulerian  convection  scheme  (At  =  0.01  in  the  calculations), 
(iii)  calculate  the  area  between  the  free  surface  and  the  x  axis,  and  (iv) 
repeat  the  above  steps.  The  procedure  described  above  is  relatively 
simple  and  does  not  require  excessive  computer  time  (about  30  min- 
utes on  IBM  3033).  However,  the  free  surface  eventually  does  become 
unstable.  The  instability  shows  up  as  a  kink  in  the  free  surface  caused 
by  the  propensity  of  the  neighboring  vortices  to  orbit  about  each  other. 
Some  of  the  instabilities  can  be  alleviated  without  smoothing.  Initial 
calculations  have  shown  that  the  vortices  near  the  y  axis  begin  to  slide 
sideways  as  the  free  surface  (or  the  vortex  sheet)  stretches.  Conse- 
quently, the  thinly  populated  regions  of  the  sheet  do  not  yield  a 
smooth  surface  and  the  flow  begins  to  leak  between  the  vortices.  Also, 
the  regions  where  the  free  surface  is  depressed  (where  the  scars 
develop)  become  overpopulated  with  vortices,  leading  to  the  growth  of 
short -wavelength  Helmholtz  instability.  This  problem  can  easily  be 
alleviated  either  by  packing  the  vortices  more  densely  near  the  y  axis 
at  the  start  of  the  calculations,  so  that  the  entire  surface  becomes 
more  or  less  uniformly  represented  at  later  times,  or  by  rediscretizing 
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the  vortex  sheet  at  suitable  time  intervals.  This  latter  technique  has 
been  used  in  the  present  calculations  (for  additional  details  see 
Sarpkaya  and  Shoaff  [Ref.  21]). 

B.     RESULTS  OF  THE  NUMERICAL  CALCULATIONS 

Numerous  runs  were  made  at  various  Froude  numbers  from  0.35 
to  0.65.  In  each  case,  the  calculations  were  carried  out  until  the  insta- 
bilities were  strong  enough  to  cause  a  numerical  floating  point  over- 
flow error  at  some  point  during  the  execution  of  the  program.  A 
numerical  filter  was  not  used  to  eliminate  or  to  delay  the  instabilities. 
Nondimensional  time  was  used  as  well  as  nondimensional  lengths  and 
speeds.  For  the  sample  results  shown  here,  the  counter-rotating  vor- 
tices were  placed  at  d0/b0  =  3.  The  start  of  their  motion,  measured  by 
T  =  (V0t/b0  -  d0 /b0),  was  T  =  -3.  The  free  surface  was  not  allowed  to 
start  to  deform  until  the  vortex  pair  was  within  2.5  units  of  the  free 
surface.  Horizontal  placement  of  the  vortex  pair  was  at  -0.5  and  0.5 
units  so  that  the  distance  between  the  two  would  be  b0. 

Figures  5  through  8  show  the  evolution  of  the  free  surface  at  rep- 
resentative times  for  a  Froude  number  of  0.53.  The  corresponding 
velocity- vector  plots  are  shown  in  Figures  9  through  12.  The  path  of 
the  main  vortex  follows  the  path  of  a  vortex  pair  approaching  a  rigid 
boundary  until  about  T  =  -1.6.  Then  it  starts  to  diverge  and,  instead  of 
turning  to  the  right  as  in  the  case  of  the  rigid  wall,  it  starts  to  move  to 
the  left  very  slowly.  This  is  because  the  strengths  of  the  discrete  vor- 
tices placed  near  the  center  line  are  inducing  a  horizontal  velocity 
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Figure  5.  Evolution  of  Free  Surface  at  T  =  0.95 
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Figure  6.  Evolution  of  Free  Surface  at  T  =  0.85 
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Figure  7.  Evolution  of  Free  Surface  at  T  ■  0.75 
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Figure  8.  Evolution  of  Free  Surface  at  T  =  0.65 
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Figure  9.  Velocity  Field  at  T  =  0.95 


25 


.lX    .)    A    J    U    ^    s\    si    n!    \l    \1    \l    \1    n! 


o 

DO 


I 


_,! 


.I  A©i  ^     J        l 


tn 


m 


^f 


n 


0 


>-    A 
<  <: 


fc  *  ^  ^  ^ 

^  V  y  x  < 

^  V  V  r  -c 

N  h  V  "V  x-  "C  < 

N  h  V  V  v  <c  < 

A  N  V*  V  v  v  "^ 

A  N  V  V  V  v  v 

^  ^  ^  v  v  v  v 

A  N  N  X-  V  \-  T 

h  ^  s  v-  v  v  v 

a  n  r^  v-  v  v  v 

a  n  is  v-  \^  x-  v 

A  N  N  S  X-  v  X- 

a  n  k  t^  ^  v  x- 

A  N  N  N  V-  ^  v 

A  h  h  h  V  \-  X- 

A  M  |\  fs  ^  t^  M 


-1 

V 
V 

A    A 
A    A 

/-  /-  A-  A 

Z.    Z.    Z-  A 

<  -d    ^  ^ 

<  <    <£  ^ 

^  -c  <  < 

V  <c   -c  < 

T    X-    «C  < 

X    V    V  < 

V  V    V  ^ 

\-  v  v  v 

X-    V   V  v 

V  V    V  T 
X-    V    V  V 


a  A 

A    A 


z.  Z-  A  A 

/L   £-   4-   I- 
£.£.£-£. 

<  Z\    Z.    Z. 
<.    -<C    <L    Z. 

<  <  <  *c 

<  <  <  < 

-c  *c  <  < 

■st  r:  <  < 

v  v  -c  < 

V    V    T    «c 


V 
V 
V 
V 
V 

A 
A 
A 
A 
A 
A 

Z. 
ZL 

Z.  Z. 

<£.  ZL 

<  <L 

<  < 

<  < 

<  < 


A     k 

A  A 

A  A 

A  A 

A  A 

Z-  A 

Z.  Z. 

^  ^L 

<£  Z^ 

<i  «£ 

<  «C 

<  <. 


A  A 

Z-  Z- 

Z.  Z- 

Z-  Z- 

Z.  Z- 


2  3 

X/BO 


T 
4 


A  A 

A  A 

-^  A 

^  /- 

Z-  /- 

/-  Z- 

/L  Z. 


Z.    Z- 
/-    Z- 


T 


5 


Figure  10.  Velocity  Field  at  T  =  0.85 
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Figure  11.  Velocity  Field  at  T  =  0.75 
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Figure  12.  Velocity  Field  at  T  =  0.65 
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component  to  the  left.  This  eventually  would  lead  to  the  vortex  moving 
up  into  the  domed  area  underneath  the  free  surface  if  the  calculations 
had  proceeded  further.  Figure  13  shows  a  sample  plot  of  the  stream- 
lines, corresponding  to  the  maximum  time  at  which  the  calculations 
remained  stable.  It  is  important  to  note  that  the  results  cited  above 
have  not  been  subjected  to  any  kind  of  numerical  filtering.  The  cal- 
culations are  generally  terminated  when  the  free  surface  forms  a  sharp 
corner  or  when  difficulties  are  encountered  in  obtaining  rapid  conver- 
gence. There  does  not  seem  to  be  any  direct  relation  between  the 
breakdown  of  the  numerical  model  and  the  attainment  of  maximum 
scar  conditions.  This  is  partly  because  the  effect  of  viscous  and  turbu- 
lent diffusion  on  the  counter-rotating  vortices  and  on  the  scars  and 
striations  is  not  yet  taken  into  account.  Nevertheless,  the  idealized 
model  gives  some  indication  of  the  surface  deformation  and  validates 
the  conclusion  that  the  rise  of  the  vortex  pair  does  not  give  rise  to 
waves  on  either  side  of  the  scars,  as  anticipated  by  Kochin,  et  al. 
[Ref.  15].  Finally,  it  is  clear  that  a  two-dimensional  numerical  model 
cannot  predict  the  existence  of  striations  which  are  essentially  three- 
dimensional  instabilities.  It  is  quite  possible  that  the  most  important 
physical  phenomenon  leading  to  the  existence  and  explanation  of  the 
dark  narrow  regions  of  wake  behind  the  ships  may  be  the  striations  or 
the  interaction  of  striations  with  scars. 
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Figure  13.  Streamlines  at  T  =  0.65 
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III.  EXPERIMENTS  AND  COMPARISONS 

A.   EXPERIMENTAL  APPARATUS 

Experiments  were  conducted  in  a  water  basin  which  was  twelve 
feet  long,  three  feet  wide,  and  four  feet  deep.  Additional  details  of  the 
basin  are  described  by  Leeker  [Ref.  22]  and  will  not  be  repeated  here. 

Two-dimensional  trailing  vortices  were  generated  through  the  use 
of  three  different  methods: 

1.  The  use  of  two  counter-rotating  plates  which,  when  fully  rotated, 
were  flush  with  the  bottom  of  the  tank,  i.e.,  they  have  "vanished" 
at  the  end  of  their  rotation,  as  far  as  the  resulting  fluid  motion  is 
concerned  (see  Figure  14),  and  no  secondary  vortices  were 
generated. 

2.  The  use  of  a  two-dimensional,  sharp-edged  slot  on  the  wall  of  a 
chamber,  attached  to  the  bottom  of  the  basin.  The  fluid  in  the 
chamber  is  driven  by  a  piston,  attached  to  a  pulley  and  weight 
system. 

3.  The  use  of  a  streamlined  nozzle  attached  to  the  aforementioned 
chamber,  directly  in  front  of  the  nozzle.  Most  of  the  experiments 
were  conducted  through  the  use  of  the  nozzle  because  of  the 
superior  quality  of  the  resulting  vortices. 

To  visualize  the  flow  created  by  the  vortices,  a  fluorescein  dye  was 
introduced  either  from  the  false  bottom  at  about  the  midpoint  of  each 
rotating  plate  or  from  the  inner  walls  of  the  nozzle.  The  dye  was  con- 
trolled by  making  it  slightly  buoyant  with  the  use  of  alcohol  and  by 
using  a  flow  control  valve  just  outside  the  tank  near  the  dye  exit. 
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Figure  14.  Rotating  Plates  and  the  Kelvin  Oval 
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The  data  collection  system  consisted  of  a  video  camera,  a  VHS 
video  cassette  recorder,  a  monitor,  and  a  date/timer.  The  date/timer 
was  used  to  imprint  on  the  screen  and  the  tape  a  running  clock  and 
date  (which  was  used  as  a  run  number)  for  each  experiment. 

B.  PROCEDURES 

The  tank  was  first  filled  to  a  suitable  level  and  allowed  to  sit  until 
the  fluid  had  come  to  rest.  Water  levels  were  varied  throughout  the 
experimentation  to  assess  the  effect  of  the  water  level  on  how  soon 
the  vortices  became  turbulent.  While  waiting  for  the  fluid  to  come  to 
rest,  the  plates  were  opened  and  the  weights  set.  After  the  fluid  was  at 
rest,  the  dye  was  slowly  introduced  into  the  system.  When  the  dye  had 
risen  just  above  the  tips  of  the  plates,  the  collimated  light  source  was 
turned  on  and  the  video  recording  system  was  set  to  run.  The  weights 
were  then  actuated,  causing  the  plates  to  rotate,  which  created  the 
two-dimensional  vortices.  The  vortices  rose  by  their  mutual  induction 
velocity.  As  the  vortices  continued  to  rise,  the  free  surface  started  to 
deform,  leading  to  the  formation  of  scars  and  striations. 

C.  DATA  ACQUISITION  SYSTEM 

A  NEC  TI-22P  CCD  camera  was  used  to  record  the  experimental 
runs  onto  a  Panasonic  AG-6300  Video  Cassette  Recorder.  The  VHS 
recorder  is  required  because  for  every  frame  of  video  information  dis- 
played   on    the    monitor,    the   VHS    tape    stores    another   frame    of 
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information  digitally  which  can  later  be  extracted  with  the  proper 
equipment.  A  Sun1  computer  with  a  motion  analysis  VP110  system  was 
used  to  analyze  the  data.  The  motion  analysis  system  allows  the  user  to 
superimpose  on  the  video  monitor  an  outline  of  areas  that  have  the 
same  intensity  or  brightness  levels.  The  brightness  or  area  to  be  out- 
lined is  controlled  by  the  user  as  he  sets  a  threshold.  The  threshold 
can  be  varied  from  zero  (black)  to  1999  (white).  The  area  outlined 
shows  up  on  the  monitor,  allowing  the  user  fine  control  of  what  is  to 
be  picked  up  by  the  computer. 

The  outlines  are  then  turned  into  digital  data  and  retained  in  a  file 
by  pixel  location  of  the  monitor  screen.  Once  the  data  is  stored,  it  can 
be  edited  and  manipulated  to  eliminate  unwanted  images  (see  Figures 
15a  and  15b).  Further  manipulation  then  allows  the  user  to  find  cen- 
troids  of  areas,  paths  of  centroids,  and  velocity  and  acceleration  of  the 
particles  on  these  paths.  This  is  all  done  through  the  use  of  the  soft- 
ware provided  in  a  package  known  as  ExpertVision.2 

Figures  16  and  17  show  the  evolution  of  the  free  surface,  as 
determined  through  the  use  of  the  ExpertVision,  for  the  Froude  num- 
bers of  Fr  =  0.39  and  0.53.  Also  noted  in  these  figures  are  the  posi- 
tions of  the  counter-rotating  vortices.  These  will  later  be  compared 
with  those  predicted  numerically. 


*Sun  is  a  registered  trademark  of  Sun  Microsystems,  Inc. 

2ExpertVision   is    a   registered    trademark   of  Motion   Analysis 
Corporation. 
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Figure  15a.  Identification  of  the  Free  Surface  and  the  Kelvin  Oval 

(early   stages) 
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Figure  15b.  Identification  of  the  Free  Surface  and  the  Kelvin  Oval 

(later   stages) 
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Figure  16.  Evolution  of  the  Free  Surface  During  the  Rise  of 

Vortices  for  Fr  =  0.39 
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Figure  17.  Evolution  of  the  Free  Surface  During  the  Rise  of 

Vortices  for  Fr  =  0.53 
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D.      COMPARISON   OF  ANALYSIS  AND  EXPERIMENTS 

In  general,  the  numerical  model  predicted  closely  several  features 
that  occurred  in  the  experiments:  the  free  surface  shape,  particularly 
the  scar  region  and  the  smooth  shape  of  the  hump  or  domed  area  in 
the  center.  Furthermore,  the  model  predicted  that  little  would  hap- 
pen beyond  y/bQ  =  3.  This  was  found  to  be  true  for  the  stages  of  defor- 
mation up  until  the  time  when  the  vortices  lost  nearly  all  of  their 
vertical  movement.  The  model  also  predicted  the  location  of  the  scars 
to  a  high  degree  of  accuracy  as  well  as  the  slight  inward  movement  of 
the  path  of  the  vortices  during  their  vertical  ascent  as  they  are  drawn 
into  the  humped  region. 

Figures  18  through  20  show  a  comparison  of  the  variation  of  the 
local  maximum  of  the  free-surface  deformation  as  a  function  of  the 
instantaneous  position  of  the  counter  rotating  vortex  pair  for  Fr  = 
0.39,  0.42,  and  0.53,  respectively.  Where  applicable,  the  agreement 
between  the  experimental  data  and  numerical  predictions  is  quite 
satisfactory,  particularly  when  one  considers  the  fact  that  the  surface 
deformations  take  place  in  a  time  period  less  than  0.3  seconds  and  the 
fact  that  the  numerical  model  is  an  inviscid  flow  idealization  of  our 
conceptual  view  of  the  phenomenon.  The  smaller  the  number  of  vor- 
tices or  marker  elements,  the  larger  the  dome  height  and  slower  the 
onset  of  the  numerical  instability.  This  is  neither  surprising  nor 
entirely  unexpected.  The  stability  of  vortex  sheets  and,  in  particular, 
the  Helmholtz  instability  have  been  the  subject  of  intense  interest.  It  is 
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Figure  18.  Comparison  of  the  Measured  and  Predicted  Scar  Elevations 

for  Fr  =  0.39 
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Figure  19.  Comparison  of  the  Measured  and  Predicted  Scar 

Elevations  for  Fr  =  0.42 
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Figure  20.  Comparison  of  the  Measured  and  Predicted  Scar 

Elevations  for  Fr  =  0.53 
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a  well-known  fact  that  an  infinitesimal  disturbance  of  wavelength  X  on 
a  plane  sheet  of  strength  y  grows  like  expfaytA),  according  to  which 
the  shorter  waves  grow  faster  [Refs.  18,  23-25].  Clearly,  the  larger  the 
number  of  vortices,  the  smaller  the  wave  length  and  the  faster  and 
sooner  the  inception  of  instability,  as  demonstrated  once  again  by  the 
present  calculations.  As  noted  earlier,  the  instability  could  have  been 
further  delayed  by  filtering  the  numerical  errors  leading  to  the  saw- 
tooth instability.  In  fact,  this  could  have  yielded  scar  elevations  even 
larger  than  those  observed  experimentally.  But  the  use  of  an  artificial 
filter  has  been  disregarded.  The  purpose  of  the  present  calculations 
was  not  the  exact  prediction  of  the  ultimate  scar  shape  but  rather  a 
search  for  new  insights  and  inspiration  in  the  simulation  of  a  complex 
three-dimensional  free-surface  flow.  The  challenge  was  to  predict  the 
rate  of  growth  in  the  face  of  ever-accumulating  errors  within  a  time 
scale  appropriate  for  the  investigation.  The  numerical  model  cannot 
and  is  not  expected  to  deal  with  three-dimensional  instabilities  and 
turbulence  seen  during  the  actual  experiments.  Although  the  model 
has  the  capability  to  vary  the  vortex  strength  with  time  to  account  for 
the  effects  of  laminar  and  turbulent  diffusion,  this  was  not  done.  By 
letting  the  vortex  strength  diminish  with  time,  the  program  may  have 
run  longer  and  the  actual  height  of  the  humped  region  may  have 
agreed  more  closely  with  the  heights  seen  during  the  experiments. 
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IV.  CONCLUSIONS 

The    investigation    reported    herein    warranted    the    following 
conclusions: 

1.  The  numerical  model,  based  on  the  use  of  a  vortex  sheet,  redis- 
cretization  of  the  sheet  segments,  and  the  desingularization  of 
the  complex  velocity,  without  the  use  of  an  additional  numerical 
filter,  is  capable  of  predicting  the  rate  of  growth  of  the  free -sur- 
face shapes  and  the  location  of  the  early  stages  of  the  free  surface 
deformations.  The  ultimate  heights  of  the  scars  and  striations  are 
dictated  by  laminar  and  turbulent  diffusion  of  vortices,  surface 
impurities,  and,  to  a  lesser  extent,  by  the  initial  shape  of  the 
counter-rotating  vortex  pair.  The  calculations  are  generally  ter- 
minated when  the  free  surface  forms  a  sharp  corner  or  when  dif- 
ficulties are  encountered  in  obtaining  rapid  convergence.  There 
does  not  seem  to  be  any  direct  relation  between  the  conditions 
prevailing  at  the  time  of  breakdown  of  the  numerical  model  and 
the  physical  conditions  existing  at  the  time  of  attainment  of 
maximum  free-surface  deformation.  Clearly,  there  are  other  phe- 
nomena, such  as  turbulence  and  three-dimensional  effects,  with 
which  the  present  numerical  model  cannot  deal. 

2.  For  any  of  the  Froude  numbers  encountered  in  the  analysis  or 
experiments,  no  wave  train  was  observed  on  either  side  of  the 
scars.  This  conclusion  was  somewhat  anticipated  on  the  basis  of 
previous  calculations.  For  very  small  Froude  numbers  (say, 
Fr<0.15),  the  scars  were  very  small,  the  vortices  followed  the 
rigid -wall  path,  and  the  scars  were  slaved  to  the  vortices. 

3.  Preliminary  experiments  have  revealed  the  existence  of  striations 
even  in  an  otherwise  two-dimensional  vortex  flow.  In  other 
words,  the  striations  observed  in  earlier  investigations  at  the 
Naval  Postgraduate  School  with  an  inclined  trailing  vortex  pair 
[Ref.  26]  were  not  due  to  the  three-dimensionality  of  the  flow.  It 
appears  that  the  striations  are  a  consequence  of  gravitational  and 
centrifugal  instability,  in  both  two-  and  three-dimensional  vortical 
flows  occurring  in  the  vicinity  of  a  free  surface.  It  also  appears 
that  the  turbulent  wake-like  features  observed  in  SAR  pictures 
are  related  to  a  vortex  mechanism:  formation  of  scars  and  stria- 
tions, intersection  of  the  striation  waves  with  the  edges  of  the 
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scars,  and  the  formation  of  small-scale  vortical  patterns.  These 
features  will  be  the  subject  of  continuing  investigation. 
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